Dynamics of Biased Domain Walls and the Devaluation Mechanism 
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We study the evolution of biased domain walls in the early universe. We explicitly discuss the 
roles played by the surface tension and volume pressure in the evolution of the walls, and quantify 
their effects by looking at the collapse of spherical wall solutions. We then apply our results to a 
particular mechanism, known as the devaluation scenario, in which the dynamics of biased domain 
walls was suggested as a possible solution to the cosmological constant problem. Our results indicate 
that devaluation will in general lead to values of the cosmological constant that differ by several 
orders of magnitude from the observationally inferred value, ~ 10~^ eV. We also argue that 
the reasons behind this are not specific to a particular realization, and are expected to persist in 
any scenario of this kind, except if a low energy cut-off on the spectra of vacuum energy densities, 
of the order of the critical density at the present time, is postulated. This implies that any such 
scenario will require a fine-tuning similar to the usual one. 



I. INTRODUCTION 

Topological defects necessarily form at cosmological 
phase transitions P, Work on defects tends to fo- 
cus on cosmic strings, which are usually benign in terms 
of their cosmological consequences. Domain walls are 
comparatively neglected since they are subject to tight 
constraints, in particular the so called Zel'dovich bound 
3. Nevertheless, they can play interesting cosmological 
roles. For example, it was thought that a particular type 
of domain wall network that would be frozen by the cos- 
mological expansion could be a candidate for the dark 
energy that is thought to be accelerating the recent uni- 
verse although detailed studies have now ruled out 
this scenario 0, [g| . Another example is that of biased 
domain wall networks 0, which may provide a way 
to evade the Zel'dovich bound and are also behind the 
recently proposed devaluation scenario 0]. 

In order to analyse this and other paradigms one needs 
to understand the dynamics of domain wall networks in 
some detail. Since they are intrinsically non-linear ob- 
jects this usually requires a careful combination of ana- 
lytic modeling and high-performance numerical simula- 
tion. On the other hand, this toolkit can be ap plied not 
only to the simplest wall networks [l^, [HI [12] but also 
to more complicated ones, such as walls with junctions 
@i S ■ present report we will take a detailed 

look at the dynamics of biased domain wall networks, 
and then use our results to carry out a thorough analysis 
of the devaluation scenario. Another recent analysis of 
this type of walls can be found in Throughout this 
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paper we shall use natural units with c = h — ks = 1- 



II. 



BASIC DOMAIN WALL EVOLUTION 



The space-time history of an infinitely thin domain 
wall can be represented as 



a = 0,1,2. 



(1) 



The action describing the dynamics of such domain walls 
is given by 



(2) 



where a is the domain wall tension, (3 — det{Pab) and 
Pab — 9tJ.vX^a'^\ world-sheet metric. The equation 

of motion of the domain walls is 

(3-''^da{^PP''''x^,) + r(;„/3'^^x:;x", = . (3) 

Although this is not easy to solve in general, we may sim- 
plify the problem considerably by considering the evolu- 
tion of a spherical domain wall in a flat FRW universe 
with line element 



(4) 
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where a is the scale factor, t is the physical time and x are 
conformal spatial coordinates. The space-time trajectory 
of a spherically symmetric domain wall can be written as 

x{t,9)—q{t) (cos 6* cos 0, cos 6* sin 0, sin 0) . (5) 

with < 4> < 27r, < 9 < n independent of t. 
Spherical domain walls in Minkowski space-time (a = 1) 
have a conserved wall invariant area [lH given by S" = 
A-n-jq^ (with the usual definition of the Lorentz factor, 
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(1 — u^) with V = adq/dt) which implies that 



R = 7^^^|<7| is a conserved waU invariant radius (this is 
analogous to Rg — "flql for a cosmic string loop). It is 
then straightforward to show that the spherical domain 
wall equation of motion in Minkowski space is given by 



— = 1 - — 



(6) 



For the more generic case with Hubble expansion the 
wall invariant radius includes the scale factor, R — 
^^^"^aq^ and the effect of the Hubble damping term follows 
trivially from the fact that the domain wall momentum 
per unit comoving area for an infinite planar domain wall 
is proportional to (so that oc a~^), leading to 



dv 



(1 



R 



iHv 



(7) 



where, H — {da/dt)/a is the Hubble parameter and we 
have defined f{v) — 2^^/^sign{—q). The other effect of 
expansion is that R is no longer a constant but becomes 
a dynamical quantity. It is also straightforward to show, 
for example using Eqn. that its evolution equation 
has the form 



dR f 3 , , 
— ^ [1- -v^ \ HR. 
dt \ 2 



(8) 



It is instructive to compare these with the (averaged) 
evolution equations for a whole domain wall network [12l | . 
These are written in terms of a characteristic length scale 
defined as 



(9) 



and a root-mean squared velocity v, and have the form 
dL 



dt 



dv 
'dt 



3Hv 



(10) 



(11) 



For the moment we are ignoring the network's energy 
losses. The matching is closer than might naively appear 
if one recalls that in the context of a simple one-scale 
model the correlation length is the only length scale in 
the problem and is therefore identified with the curvature 
radius [lB| , so even though the velocities in both cases are 
defined differently their evolution is the same except for 
the behavior of the phenomenological curvature param- 
eter k (which in the context of analytic modeling has a 
non-trivial velocity dependence 'W\ ) . 

On the other hand, when contrasting the length scales 
the key difference is that L is to be understood as av- 
eraged over a network (hence a sizable volume) whereas 
R in our above discussion essentially applies to a single 
defect. Let us consider a domain wall that separates do- 
mains 1 (inside) and 2 (outside), and assume that the 



volume fraction occupied by each vacuum is /i and /2 
with /i < /2. Then a definition of characteristic scale 
which is more directly related to the typical size of the 
domain walls is given by 



fi 



(12) 



In the limit /i ^ /2 we expect to have only isolated 
domains so that /1//2 oc L^g/a^ from which one obtains 



(13) 



Because of this difference it is not possible to find a simple 
expression that smoothly interpolates between the two 
limits and is valid throughout: a wall network in such an 
intermediate regime will have highly non-trivial statisti- 
cal properties. (Note that this is not a serious problem 
at all: the dynamics of biased walls is sufficiently fast to 
ensure that no network will ever be in such an intermedi- 
ate state for a significant amount of time.) Substituting 
the above relation into Eqn. [TO] we get the suggestive 



dLcff 
dt 



l--vn HL 



^eff . 



(14) 



All the above discussion in essentially classical, in partic- 
ular by neglecting energy losses through radiative mecha- 
nisms. These can still be added to the model phenomeno- 
logically when needed. However, for our present purposes 
we shall be more interested in modeling the effects of bi- 
asing, as described below. 



III. BIASED DOMAIN WALLS 

In the simplest model of biased domain walls there is 
an asymmetry between the two minima of the potential 
0, [1] . The volume pressure from the biasing provides a 
further mechanism which affects the dynamics of these 
walls. Depending on its importance relative to other 
processes — most notably the surface tension — the walls 
may be long-lived (as in the standard case) or disappear 
almost immediately. 



A. Qualitative analysis 



Let us consider a simple model 



with the tilted potential 
^(0) = ^ry^ 



— - 1 I + K- 



(15) 



(16) 
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It then follows that the height of the potential barrier 
and surface tension are respectively 



^0 = iv' 



(17) 



(18) 



while the wall thickness and the asymmetry parameter 
(or energy difference between the two vacua) are respec- 
tively 



e = 2k.Vo . 



(19) 



(20) 



Now, the surface pressure (from the tension force) and 
the volume pressure (from the energy difference between 
vacua) are 



Pv 



(21) 



(22) 



At early times the surface tension tends to dominate (due 
to the small curvature radii), but when the domains be- 
come large enough they will decay. This happens when 



e 



(23) 



For example, we phenomenologically expect that at the 
Ginzburg temperature (when the network first becomes 
well-defined) its correlation length should approximately 
be given hy L ^ (Ary)"^ and in this case we get 



Pv\ 
PtJto 



2VA 



(24) 



typically one might expect this to be smaller but not 
much less than unity, but in principle there is enough pa- 
rameter freedom to make it much larger or much smaller. 
When the volume pressure dominates, one expects that 
the walls will move with an acceleration 



(7 



On the other hand, if the surface pressure dominates ini- 
tially the walls may survive long enough to reach a linear 
scaling regime, L ~ i. How fast the volume pressure 
becomes important, and hence how fast the walls can 
decay is a key issue in several cosmological scenarios, in- 
cluding devaluation. Moreover, the Zel'dovich bound 0] 
provides an additional (and often limiting) observational 
constraint. 



B. Analytic modeling 

Modeling the effect of this bias is fairly straightfor- 
ward in the context of the analytic model we introduced 
above. Suppose that the two vacua on either side of the 
domain wall have different energy densities VSn and Vout 
and define e — — Vout ■ The resulting effect on domain 
wall dynamics is very easy to understand if we consider 
a planar domain wall in Minkowski space. In this case 
the inner and outer regions of the domain wall are not 
well defined but we shall assume that the domain wall is 
moving in the direction out —>■ in. Energy conservation 
implies that 



d((T7) = vdte , 



from which we get 



dv 
'dt 



(26) 



(27) 



we can immediately note that this is coincides with Eqn. 
(PS)) apart from the relativistic gamma factors. With our 
conventions, if e > then the wall gains momentum. We 
thus confirm the expectation that the domain wall will 
feel a pressure which will tend to drive it into the region 
of higher energy density. We can alternatively write this 
as 



dv 



„2^3/2 



-f] 



Rv 



(28) 



which makes it clear that e acts like an effective cur- 
vature, which therefore accelerates the wall. Its distin- 
guishing feature is that the scale is set by the micro- 
physics of the model in question (ultimately by the form 
of the potential), rather than the macroscopic dynamics 
as is the case for the usual curvature radius of the walls. 
Moreover, this length scale is constant, whereas the usual 
curvature radius increases as the walls evolve. This im- 
plies that the volume pressure term gradually becomes 
more important. Having said that, note that the extra 7 
factor can switch this term off if the walls become ultra- 
relativistic {v ^ 1). This discussion makes it clear that 
we can now add this volume pressure correction to our 
domain wall evolution equation ([7]), yielding 



(25) '3'" equivalently 



dv 
'dt 



(1 



R 



m 

R 



2,Hv 

cr7 



7 
Rv 



SHv 



(29) 



(30) 



Taking into account the effect of the bias in Eqn. (O, 
which describes the evolution of the wall invariant radius, 
i?, we obtain 

^=(i-IvAhR+'^, (31) 
dt \ 2 ) 2(77 ' 
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or again using the definition of Ry provided in Eqn. [28] 
it can be written in the more suggestive form 

which should be compared with Eqn. Notice that 

this additional term due to the bias is associated with an 
extra energy loss by the network. As discussed above, the 
importance of this term is expected to grow as the wall 
evolves, although it switches off in the ultra-relativistic 
limit. 

It obviously follows from the above discussion that 
this model reproduces the result Q that domains with 
a larger energy density will decay when their typical size 
R > cr/e, and indeed it provides a more quantitative 
estimate — note that the relativistic gamma factor may 
be significant. We may also compare the importance 
of the pressure term with that of the Hubble damping 
term in Eqn. ([7]). For non-relativistic domain walls the 
pressure term dominates over the Hubble damping term 
slightly earlier than over the curvature term (assuming 
that R < However, for R and v not too 

small the two criteria are very similar. If i? < and 
the above criteria are satisfied then the domains with a 
larger energy density disappear exponentially fast. This 
result is the basis of the devaluation mechanism. 



C. A numerical example 

As an illustrative example of the effect of a bias on 
the dynamics of the domain walls, we have solved Eqns. 
([29H3T|) numerically in the radiation era for a spherical 
domain wall, initially at rest, with a larger vacuum en- 
ergy density in its inner domain. This allows us to deter- 
mine the evolution of the domain wall's physical radius, 
aq = i?7~^/^, until collapse (which was defined as the 
moment when q vanishes). In this section we shall as- 
sume that Gi — a{ti) ~ 1 where ti = H^^ /2 is the initial 
time (we take ti — 1). 

Fig. [T] shows the time of collapse of a spherical domain 
wall and the relative difference between these collapse 
times as a function of the ratio between the initial wall 
invariant radius, i?, and the initial Hubble radius. The 
bias term acts as a further mechanism which accelerates 
the domain walls, allowing them to overcome the Hubble 
damping term faster, and hence making them collapse in 
a shorter period of time. The relative importance of this 
effect grows as the initial radius increases. However, the 
curvature term in the equations becomes negligible for 
large R and therefore in this limit the collapse of the wall 
is essentially determined by the bias term. As a result, 
in this limit the relative difference tends to a constant. 

Fig. [2] shows the relative difference between the col- 
lapse time, with and without an energy difference be- 
tween minima, for e/tr ranging from to 100 (for a spher- 
ical domain wall with Ri = H~ ^ ) . Notice that as the bias 
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FIG. 1: The top panel shows the collapse time of a spherical 
domain wall, with and without an energy difference between 
the inner and outer domains, as function of the ratio between 
the initial wall invariant radius and the initial horizon. Time 
is in units of initial time, ti. The lower panel depicts the rela- 
tive difference between the time of collapse with and without 
bias. In both cases we have assumed e/cr = 1. 
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FIG. 2: The relative difference between time of collapse with 
and without bias, for a spherical domain wall with Ri — H^^ , 
for e/a ranging from to 100. 



term's importance grows, this relative difference tends to 
a constant. Since the bias term makes the velocity ultra- 
relativistic faster, the time of collapse with bias has an 
obvious lower bound, given by the time a photon takes to 
travel a distance equal to the wall's physical radius. In a 
flat FRW universe: tj?"™ = (1 + 9^/2)^ For this particu- 
lar case, this yields t™™ ~ 4, and the relative difference 
approaches A ~ 0.33 as the ratio e/a increases. 



IV. THE DEVALUATION MECHANISM 

Devaluation Q is a mechanism which aims to explain 
the observed small non-zero value of the cosmological 
constant using the dynamics of biased domain walls. The 
key physical idea is that under plausible circumstances a 
network of unstable domain walls might form at a certain 
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critical temperature, dividing the universe in many dif- 
ferent regions with different values of the vacuum energy 
density. Then domain walls separating different vacuum 
domains will feel a pressure which will tend to suppress 
those with higher energy thus driving the universe to- 
wards lower and lower values of the vacuum energy den- 
sity. 



A. A simple realization 



A simple potential that has the relevant features is 

F(</>) = Focos(^^) -e^-HC, (33) 

where Vb is effectively the barrier height and e is the en- 
ergy difference between minima. From the original paper 
0, the expectation is that 



V'^Vq^'^ > lOMeV 



and 



.1/4 



plii^ 10-' eY. 



(34) 



(35) 



It is important to realize that, in this particular model, 
the domain walls interpolating between successive pairs 
of minima in the above potential all have similar tensions. 
The model also requires a large number of minima, but 
the exact number is actually not relevant for the analy- 
sis that follows. It was argued that domains with neg- 
ative energy vacua would be suppressed so that domain 
wall dynamics would lower the cosmological constant to a 
value just below zero {p^ac ^ e ^ Pcrit)- We will question 
the validity of this assumption later on. 

The initial conditions are expected to be such that 
?7^i, which is significantly smaller 



-1/4 



than the Hubble radius. 



mpi 



Vr 



1/4 



< 10^ 



(36) 



Clearly, with Li ^ rj^^ , the surface pressure pT ^ rf ^ 
Vo will initially dominate the volume pressure pv ~ e and 
the walls will initially be very stable. So for this choice 
of parameters, devaluation is not very efficient. 

On the other hand, the largest correlation length we 
can have (which corresponds to the largest instability) is 
L ^ t, and recalling that 



t ^ 



in the radiation era and 



mpi 
2^2 



mpi 



2.3/2yV2 



(37) 



(38) 



in the matter era, we find that in this maximally instabil- 
ity case the domain walls become unstable and disappear 
will for 



T 
mpi 



V 



Vo mpi 



\mpi 



3/2 



77 / mpi \ 



Vq mpi \ Teq 



(39) 



(40) 



For the decay to happen in the radiation era we need 

T>Te, ^ 



L eg — 1 eV, and using Eq. ((35|l we find 
77™<j < lOOkeV, 



(41) 



which is clearly incompatible with the assumed bound 
given by Eq. p4)) . It follows that if this mechanism 
(at least as originally presented) is to explain the ob- 
served value of the cosmological constant, the domain 
walls must survive through the radiation era, and decay 
only in the matter era. In other words, in the original 
scenario devaluation must occur late in the history of the 
universe, and not early. This is due to the energy scale 
required to match the observed value of the cosmological 
constant, and is another manifestation of the underlying 
fine-tuning. Repeating the calculation for a decay during 
the matter era (and ignoring the effect of the recent accel- 
eration phase on the expansion rate, which is negligible 
in this context) we now find 



Vmat < lOMeV, 



(42) 



which saturates the bound given by Eq. p4|) . 

Therefore the best we can do is to have a network that 
disappears around today. If so, and again neglecting the 
effect of the recent dark energy domination, we would 
expect the cosmological constant to be 



1/4 
Pvac 

mpi 



( V 



\mpi 



3/4 



mpi 



3/8 



eq 



\mpi 



1/8 



and since Tf.q ~ 1 eV and Tq ^ 2 x 10 eV we get 



rvac 



> 3/4 

^ ^ 10^3 eV. 



mpi 



(43) 



(44) 



from which we would get plit — 10^"^ eV for 77 ~ 10 MeV 
as previously stated. 

There is, however, an obvious problem with such a 
scenario: 10 MeV domain walls decaying today are ob- 
servationally ruled out. For the rather classic Zel'dovich 

bound [3, ~ 1 MeV, we only get pHt ~ lO""* eV. This 
is optimistic both in the sense that the observational 
bound is somewhat lower, and that we really want the 
network to decay a bit before today to be clear of obser- 
vational problems. 

Indeed, for late devaluation there are several extra 
requirements (not present for early devaluation) which 
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need to be satisfied. If devaluation is not complete and 
there are some walls still around, their average contri- 
bution to the energy density of the universe needs to be 
Pwaiis < 10~^/9crit- Otherwise, assuming that the charac- 
teristic size of the domains is > /lOO, there would be 
a detectable contribution of domain walls to the cosmic 
microwave background anisotropies. This means that the 
wall tensions are strongly constrained. Wc know from 
Eqn. (|23|) that if we have two contiguous domains with 
a vacuum energy difference of e then the domain with 
larger vacuum density will be exponentially suppressed 
when e > cr/L. But a/L is the average energy density 
of the domain wall network which, as we saw, must be 
at least five orders of magnitude smaller than the critical 
density. So in this case, we obtain 

pV^<10-4eV, (45) 

and devaluation would lead to a vacuum energy den- 
sity significantly smaller than the critical density at the 
present time. In fact this number might even be smaller, 
depending on the domain wall tensions. 



B. General features of the mechanism 

Beyond particular realizations of the devaluation 
mechanism, several general comments can be made about 
the underlying physical scenario. First, let us note that 
this simple model may allow for positive and negative 
values of the vacuum energy density. It is argued that 
gravitation itself may prevent the vacuum energy density 
from attaining negative values. However, here we will 
be mainly concerned with the evolution of domain wall 
networks during the matter and radiation eras in which 
the contribution of the domain wall and the vacuum en- 
ergy densities can be neglected. In this context, there 
is no cutoff that prevents the vacuum energy densities 
from attaining negative values. Such a mechanism can 
only operate when the vacuum energy densities are the 
dominant contribution for the dynamics of the universe. 
Hence, a low-energy cutoff of the order of the critical 
density at the present time needs to be introduced by 
hand in the devaluation model, which is clearly not an 
attractive feature of the model. 

In this scenario a domain wall network forms at a criti- 
cal temperature Tc ^ Vq^^- It is straightforward to put a 
lower bound the number, N, of domains that are initially 
present in a region with a comoving size of the order of 
the Hubble radius at the present day. 

Assuming a fixed energy difference between successive 
minima and that the barrier height and tuning of the 
potential are roughly comparable as in the simple de- 
valuation toy model introduced above we see that the 



difference in energy densities between successive minima 
is bounded from below 

e<^-Tc(ToTcq)3/2 (47) 

if all the minima are populated at the time when the 
network forms. In fact, if we assume that close to Tc the 
domains rapidly attain a typical size of the order of the 
Hubble radius then 

e-T,(Toreq)'/' (48) 

which is much larger that the energy density at the 
present time. Of course, not all the minima need to be 
populated and in this case e can be smaller. However, in 
this case the domain wall tension is no longer the same 
for all domain walls since domain walls will in general 
interpolate between distant minima. We also note that 
a single domain with the lowest possible energy density 
does not necessarily survive domain wall evolution since, 
even in the absence of the devaluation mechanism (which 
will only operate for e > a/L), domain wall dynamics in 
the scaling regime naturally leads to the suppression of 
most of the available domains each Hubble time. 

The above analysis also confirms the naive expectation 
that the devaluation mechanism can only be effective if 
L < H^^ . If this is not the case (either because the 
walls are somehow pushed outside the horizon, or in the 
opposite limit, L <^ H~^) the dynamics of the domain 
walls is even less efficient in suppressing domains with 
larger values of the vacuum density. In the absence of 
friction, and assuming that L < -ff^^, a domain wall 
network will in fact approach a scaling regime with L oc 
H~\ which is quite generically an attractor solution [13, 
lia, hj. Friction may slow down the domain walls, and 
if they are sufRciently light the evolution may be friction 
dominated up to the present time, but this will not help. 

Therefore the devaluation scenario does not naturally 
lead to the required value of the vacuum energy density. 
A simple and physically intuitive way of expressing this 
is in terms of fine-tuning. A cosmological constant in 
the standard scenario is considered unappealing because 
the value that observations indicate for its magnitude is 
many orders of magnitude below what one would natu- 
rally expect from particle physics considerations. From 
this perspective the motivation of the devaluation sce- 
nario is that it would lead to an otherwise extremely 
small value in a natural way. However, this is not so 
because the same fine-tuning problem is still there. 

Of course many other features of the simplest imple- 
mentation of devaluation can be relaxed. We may have 
complex domain wall networks with junctions, domain 
walls with tensions that may be correlated with the dif- 
ferences between the vacuum energy densities. In fact 
the domain walls may not all be formed at the same time 
with additional domains with smaller energy densities 
separated by low tension domain walls being formed only 
at smaller critical temperatures. Still, it is clear that the 
devaluation mechanism is too efficient and consequently 
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it will always be necessary to introduce by hand a low 
energy cut-ofF of the order of the critical density at the 
present day. 

V. DISCUSSION 

In more realistic particle physics scenarios there will 
typically be many coupled scalar fields that can lead to 
domain walls. From a phenomenological point of view, 
there are several reasons why this more general case may 
differ from the simplest implementation of the devalua- 
tion mechanism. The potential may be significantly more 
complicated, as in landscape-type scenarios. This by it- 
self need not be a great advantage, since in any case 
energy minimization criteria will always favor evolution 
down the potential, regardless of the number of fields. (In 
particular, any number of uncoupled fields will lead to a 
scenario very similar to the single-field case.) More im- 
portant, though, is the existence of coupled fields, since 
this generically leads to domain wall networks with junc- 
tions, and also to more complicated spectra of wall ten- 
sions, which need not necessarily be of comparable mag- 
nitude. 

The presence of walls with significantly different ten- 
sions may be important for the network dynamics. Re- 
call that biased walls will decay when their characteristic 
size grows to L > cr/e. Since in this model e is effectively 
fixed to the observed vacuum energy density, we see that, 
according to this criterium, higher-tension walls are actu- 
ally more stable than lower-tension walls (since a larger 
characteristic size is needed to make them decay). Notice 
that this is contrary to standard energy minimization ar- 
guments, whereby higher-tension walls tend to decay into 
lower-tension walls. 

If we now consider a network of walls with junctions 
and a non-trivial hierarchy of tensions, and temporarily 
assume that the various types of walls with different ten- 



sions have comparable characteristic sizes (which may be 
too naive an assumption), then one may reach a threshold 
where the lowest-tension walls present become unstable 
and decay. This may render the whole network unsta- 
ble and make it disappear well below what one would 
expect from a stability analysis for the higher-tension 
walls. Therefore this mechanism may increase the effi- 
ciency of devaluation. However, given that the devalu- 
ation mechanism is already too efficient in its simplest 
implementation this extra efficiency will not help. Hence 
a low energy cut-off of the order of the critical density 
at the present day is needed in order for devaluation to 
stop at the observed value of the dark energy density. 

Although in this article we have dealt with biased do- 
main walls in a post-inflationary context, it is important 
to realize that in the context of string landscape scenarios 
they can arise during the inflationary phase itself, as the 
inflationary field samples multiple non-degenerate vacua 
From our point of view this would be an exam- 
ple of early devaluation, with the additional difference 
that the energy scales involved are much higher (and the 
walls would typically be horizon-sized). The subsequent 
evolution of the system will depend on several factors, 
including the relative populations of the different vacua 
and the scales of the various deviations from vacuum de- 
generacy. Requiring that the standard cosmological be- 
havior is recovered by the time of nucleosynthesis leads 
to non-trivial constraints on the model parameters. 
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